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An  image  analysis  software  system  has  been  developed  for  a  UNIX 
operating  system  and  a  Tektronix  terminal.  It  is  designed  for  the 
particles  and/or  pores  in  construction  materials.  This  system  provides 
a  means  to  measure  many  discrete  fields  of  view,  in  order  to  obtain 
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either  an  electron  or  light  microscope,  or  from  an  ordinary  camera. 
Either  negatives  or  positives  can  be  used.  A  new  method  based  on  fuzzy 
probability  has  been  developed  to  segment  digital  images  into  binary 
images.  It  is  more  consistent  than  existing  techniques.  The  accuracy 
of  the  system  has  been  verified  with  images  having  known  geometric 
properties.  Measurements  have  an  error  of  less  than  1%  with  an  image 
of  appropriate  resolution.  The  system  makes  different  measurements  on 
separate  particles  and  cut  planes  through  massive  samples.  In  the 
former  case,  the  measurements  are:  the  perimeters,  areas,  and  maximum 
chords  of  particles  in  orthogonal  directions,  the  maximum,  minimum  and 
mean  values  of  these  parameters,  and  their  distributions.  Also,  the 
area  of  the  image,  the  area  faction  of  particles  and  the  number  of 
particles  are  measured.  In  the  case  of  the  cut  planes,  the  estimated 
parameters  arei  volume  fraction  of  particles,  surface  area  per  unit 
bulk  volume,  and  surface  area  per  unit  particle  volume.  The  size  dis- 
tribution, total  number  of  particles  and  mean  diameter  can  also  be 
estimated . 
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INTRODUCTION 

The  goal  of  quantitative  microscopy  is  to  measure  the 
geometric  properties  of  particles  or  mi c ros t ru c t ure s .  These 
measurements  are  performed  on  two  dimensional  images.  Stereology 
is  used  to  estimate  the  three  dimensional  properties  from  the 
direct  measurements  made  on  the  two  dimensional  images.  Manual 
microscopic  examinations  are  time  consuming  and  tedious.  In 
order  to  obtain  statistically  meaningful  results,  it  is  usually 
required  that  many  measurements  be  made.  This  requirement  has 
been  ignored  frequently  because  of  the  effort  that  would  be 
required.  Recently,  computer  image  analysis  has  been  applied  to 
quantitative  microscopy.  Because  the  major  advantage  of  a  com- 
puter is  its  capability  for  processing  data,  computer  image 
analysis  makes  practical  the  measurement  of  many  objects  in  many 
images.  Therefore,  statistical  requirements  can  be  satisfied 
with  little  human  effort. 


The  process  of  quantitative  image  analysis  is  as  follows. 
The  original,  analog  image  is  first  transformed  into  a  digital 
image  by  a  digitizer.  A  digital  image  is  a  two  dimensional  array 
of  numbers,  each  number  representing  the  light  intensity,  or  gray 
level,  of  the  corresponding  point  in  the  original  image.  The 
points  are  usually  called  pixels.  In  order  to  display  the  image 
on  a  typical  computer  monitor,  the  digital  image  must  next  be 
segmented  into  a  binary  image.  A  binary  image  consists  of  two 
types  of  pixels:  object  pixels  and  background  pixels.  Next,  it 
is   sometimes   possible   to   improve  the  binary  image  by  removing 


imperfections  and  separateing  touching  objects.  This  improvement 
is  based  upon  the  user's  knowledge  of  what  the  real  image  looks 
like.  Finally,  the  desired  geometric  properties  of  the  objects 
of  interest  are  measured  by  the  computer  on  the  segmented,  and 
improved,  image. 

Some  images  are  of  discrete  particles.  In  this  case,  one  can 
view  the  full  profile  of  each  particle.  When  this  is  the  case, 
the  following  measurements  can  be  made  on  the  image. 

1.  Total  number  of  particles 

2.  Number  of  particles  per  unit  image  area 

3.  Projected  perimeters  of  particles 

4.  Projected  areas  of  particles 

5.  Projected  maximum  chords  of  particles 

6.  Distribution  of  projected  perimeters 

7.  Distribution  of  projected  areas 

8.  Distribution  of  projected  maximum  chords 

9.  Maximal,  minimal  and  mean  values  of  projected  perimeters 

10.  Maximal,  minimal  and  mean  values  of  projected  areas 

11.  Maximal,  minimal  and  mean  values  of  projected  maximum  chords 

Other  images  are  of  a  plane  cut  through  a  solid  mass  of 
material.  Such  a  plane  will  generally  intersect  imbedded  objects 
that  are  of  interest.  The  following  parameters  can  be  determined 
from  measurements  of  these  intersected  objects. 


1.   Volume  fraction  of  objects 


2.  Surface  area  of  objects  per  unit  bulk  volume 

3.  Surface  area  of  objects  per  unit  object  volume 

4.  Size  distribution  of  objects  s 

5.  Number  of  objects  per  unit  bulk  volume 

6.  Mean  size  of  objects 

7.  Mean  value  of  the  maximum  dimensions  of  objects 

Many  image  analyzers  are  not  specially  designed  for  the 
study  of  materials.  There  are  several  common  computer  methods 
for  segmenting  an  image  that  are  not  specially  designed  for  quan- 
titative image  analysis  because  this  requires  more  accurate  seg- 
mentation than  non-quantitative  image  analysis.  Also,  object 
labeling  and  counting  are  important  in  quantitative  microscopy. 
This  is  very  difficult  for  some  automatic  image  analysis  systems 
because  it  involves  the  process  of  individual  object  labeling. 
Sometimes  geometric  properties  such  as  the  areas  and  perimeters 
of  individual  objects  are  desired.  These  individual  measurements 
require  object  labeling  and  some  image  analyzers  cannot  do  this. 


The  current  research  work  has  developed  an  automatic  quanti- 
tative image  analysis  (AQIA)  software  system  specially  designed 
for  materials'  analysis.  It  is  implemented  on  an  UNIX  operating 
system  with  a  Tektronix  graphics  terminal  as  a  monitor.  Compared 
with  many  image  analyzers,  it  is  almost  hardware  independent. 
With  a  few  changes,  it  can  be  run  on  any  UNIX  operating  system 
with  any  high  resolution  monochromatic  monitor.  With  some  addi- 
tional changes  it  can  also  be  transported  to  a  microcomputer  with 
an  appropriate  graphic   terminal.   Therefore,   it   is   much   less 


expensive  than  the  existing  hardware-type  automatic  image 
analyzers.  It  is  specially  developed  for  quantitative  microscopy 
of  construction  materials. 

This  system  tries  to  combine  the  advantages  of  manual  exami- 
nation and  automatic  image  analysis.  It  has  some  unique  features. 
It  simulates  the  process  of  human  segmentation  to  make  this  pro- 
cess more  realistic.  It  labels  objects  to  make  individual  meas- 
urements easy.  With  labeled  objects,  AQIA  can  process  an  image 
more  accurately  and  with  more  flexibility  than  can  a  non-labeling 
analysis  system.  Because  the  reliability  of  image  measurement 
depends  on  the  accuracy  of  image  segmentation  and  processing, 
AQIA  is  expected  to  provide  more  reliable  measurements. 


AQIA  can  have  a  variety  of  specific  applications  to  cement 
and  concrete  research.  Image  analysis  may  be  the  only  reliable 
technique  for  some  measurements  such  as  the  distribution  of  big 
pores  in  concrete  and  flyash  particle  sizes.  AQIA  can  provide 
solutions  to  these  kinds  of  difficult  measurements.  However,  the 
aim  in  the  development  of  the  AQIA  system  has  been  to  produce  a 
package  of  programs  that  are  generally  applicable  to  the  broad 
spectrum  of  construction  materials  research.  Thus,  specialized 
programs  that  produce  only  one  result  have  been  avoided.  It  is 
hoped  that  this  has  produced  a  package  that  will  be  broadly 
applicable  to  researchers. 


BRIEF  GUIDE  TO  THE  AQIA  PROGRAMS 

The  AQIA  system  is  a  package  of  13  major  programs  and  6  aux- 
iliary programs.  The  following  is  a  brief  description  of  the  use 
of  these  in  a  typical  image  processing  application.  The  program 
names  are  written  in  capital  letters,  and  are  listed  in  the  order 
in  which  they  would  be  used  usually. 

First,  DSCRB  initiates  the  package.  It  creates  a  file  to 
hold  the  processed  image.  In  the  file,  the  initiation  time  will 
be  recorded  automatically.  DSCRB  also  allows  the  user  to  briefly 
describe  the  sample  from  which  the  image  has  been  taken,  and  to 
enter  the  length  of  one  pixel.  The  segmented  image  and  other 
relevant  information  will  be  appended  to  this  file  as  processing 
continues. 


REDUCE  extracts  the  image  gray  level  data  from  the  raw  image 
file  produced  by  the  digitizer.  It  also  allows  the  user  to 
reduce  the  image  size  and/or  the  number  of  gray  levels.  REDUCE 
produces  a  reduced  image.  After  REDUCE,  the  user  must  select  the 
method  of  segmentation. 

AQIA  has  two  ways  to  segment  a  digital  image:  straight  gray 
level  histograming  and  edge-detection  histogram  thresholding.  If 
straight  histograming  is  desired,  HISTO  will  be  called  next.  It 
compiles  a  histogram  of  the  occurrance  of  various  gray  levels. 
If  edge-detection  is  desired,  EGEHIS  will  be  called  instead.  It 
singles  out  the  boundaries  of  objects,  and  compiles  a  histogram 
of  the  occurrance  of  gray  level  gradients. 


Then,  the  user  must  decide  if  he/she  wants  automatic  selec- 
tion of  thresholds  for  either  histogram.  If  the  user  chooses  the 
automatic  method.  FUZZTHR  must  be  called  next.  FUZZTHR  will  sug- 
gest thresholding  bounds  based  on  fuzzy  set  theory. 

With  or  without  automatic  thresholding,  BITSHOW  must  be 
called  next.  BITSHOW  asks  the  user  to  enter  the  lower  and  upper 
bounds  for  thresholding.  It  then  displays  the  binary  image  on 
the  screen  with  pixels  between  the  bounds  being  illuminated.  It 
allows  the  user  to  change  the  bounds  until  he/she  is  satisfied 
with  the  binary  image.  If  there  are  some  spots  in  the  objects 
that  need  to  be  filled,  'bitshow'  allows  the  user  to  choose 
either  the  overall  or  individual  filling  mode,  or  the  individual 
followed  by  the  overall  mode.  It  outputs  the  spot-free  binary 
image  to  the  file  created  by  DSCRB. 

Next,  the  user  must  decide  if  he/she  needs  to  trim  the  bord- 
ers of  the  binary  image.  If  the  user  needs  this,  TRIM  will  be 
called  next.  If  there  is  more  than  one  image  of  the  same  sample, 
the  user  may  call  EQUAL  to  cut  all  the  images  to  the  same  size  in 
order  to  make  statistical  calculations  meaningful. 

After  this,  the  user  must  decide  if  he/she  needs  erosion  of 
the  binary  image  to  separate  objects  before  measurement.  If  ero- 
sion is  desired,  EROD.CALC  will  be  called  next.  It  allows  the 
user  to  choose  either  the  overall  or  individual  mode  of  erosion, 
or  the  overall  followed  by  the  individual  mode.  After  erosion, 
the   appropriate   dilation   will  be  performed.   Then,  the  desired 


measurements  will  be  performed  on  the  processed  binary  image. 
EROD.CALC  can  perform  measurements  on  projections  of  separate 
particles  or  sections  cut  through  a  solid  sample.  EROD.CALC  will 
output  the  final  results  of  the  measurements,  and  the 
eroded/dilated  image. 

If  no  erosion  is  needed,  the  binary  image  can  be  directly 
measured  by  the  following  three  programs.  SEP.CALC  measures  the 
projections  of  separated  particles.  CUT.CALC  estimates  3- 
dimensional  properties  of  mi cr os t rue t ure s  from  2-d imens i ona 1 
measurements  on  cut  planes.  ALL.CALC  does  both.  Finally,  NEWPSD 
can  be  used  to  estimate  the  size  distribution  of  particles  that 
are  imbedded  in  a  massive  sample  when  the  image  is  of  a  plane  cut 
through  i  t . 

For  the  user's  convenience,  AQIA  also  provides  6  auxiliary 
programs.   These  are: 

1.  DESIGN    3.  INF       5.  RESULTS 

2.  SHOW      4.  HARDCP    6.  PLOTPREP 

DESIGN  allows  the  user  to  find  the  number  of  required  images 
for  volume  fraction  estimation,  given  the  image  size,  the 
estimated  mean  value  of  the  area  fraction,  the  desired  estimation 
accuracy  (half-width)  and  the  desired  confidence  level. 

SHOW  simply  displays  a  binary  image  on  the  screen,  but  does 
not  allow  processing.  INF  gives  the  relevant  information  on  an 
image  such  as  the  size  of  the  binary  image,   the   size   reduction 
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factor,  the  number  of  gray  levels,  the  percentages  of  1  and  0 
pixels,  the  thresholding  bounds,  etc..  HARDCP  gives  the  user  a 
hard  copy  of  the  binary  image  if  the  user  is  working  on  a  Tek- 
tronix 4014  terminal  equipped  with  a  hardcopy  unit. 


RESULTS  either  displays  the  numeric  results  of  the  measuring 
programs  on  the  screen  or  prints  them  on  a  specified  printer. 
PLOTPREP  collects  the  results  of  measurements  from  EROD.CALC,  or 
SEP.CALC,  or  CUT.CALC,  or  ALL.CALC,  to  provide  input  for  plotting 
by  the  UNIX  Qplot  program. 


This  flow  chart  summarizes  the  use  of  most  of  these  programs. 


DSCRB  to  initiate| 

1 

Digital  image| 

1 

REDUCE  to  compress 
image  data 

HISTO  to  obtain 
gray  level 
h i  s  t o gram 


FUZZTHR  for 
selection 


BITSHOW  to  segment 

image  and  output 

binary  image 


EGEHIS  to  obtain 
object  boundaries  and 
his  togram 


SEP.CALC  for  particles 
CUT.CALC  for  cut  planes, 

ALL.CALC  for  both,  and 
NEWPSD  for  distribution 
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DEMONSTRATION  OF  SOME  OF  AQIA's  FEATURES 

Typical  Results  for  Separated  Particles 

The  entire  AQIA  package  has  more  features  than  can  be  demon- 
strated in  this  summary.  However,  some  of  the  more  important 
functions  will  be  described  here.  Figures  1  and  2  show  two  digi- 
tized images  of  some  particles  of  quartz  sand.  These  images  have 
already  been  prepared  for  measurement  by  separating  and  labeling 
the  particles  that  touch  (followed  by  dilation  of  the  objects). 
These  images  will  be  used  to  demonstrate  the  typical  results  of 
using  AQIA  to  measure  separate  particles. 


• 


Figure  1  Figure  2 

Images  of  Quartz  Sand  Particles 
These  two  images  are  called  Sample  1  and   Sample   2.    Since 
these   two   images  are  from  the  same  material,  the  mean  values  of 
the  measurements  from  separate  images  are   also   calculated.    In 

the   following   tables,   the  notation  "x  en"  means  that  the  value 

4 
"x"  is  multiplied  by  10  ;  the  notation,  "x  e  "  means  that  the  value 

x  is  multiplied  by  the  factor,  10  . 
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Table  1  shows  the  basic  individual  measurements  that  include 
perimeters,  areas,  and  maximum  chords  in  orthogonal  directions. 
("Extent"  means  maximum  chord).  Note  that  dabel  1  has  been 
assigned  to  the  background  so  that  there  is  no  object  labeled  1 
listed  in  the  table.  Table  1  also  gives  the  size  (area)  of  the 
whole  image,  the  area  fraction  of  objects  (intersections),  the 
total  number  of  objects,  the  number  of  objects  per  unit  image 
area,  and  the  maximum,  minimum  and  mean  values  of  the  individual 
measurements.  Sample  1  was  processed  first.  Therefore,  the  mean 
values  given  in  parentheses  are  the  same  as  the  Sample  1  values. 
When  Sample  2  was  processed  next,  the  values  in  the  parentheses 
were  updated  to  reflect  the  means  of  both  samples. 


Table  2  shows  the  distribution  of  the  perimeters.  The 
number  frequency,  the  relative  frequency  and  the  cumulative  fre- 
quency are  included.  The  interval  between  classes  is  300 
microns.  Taking  the  data  for  the  4th  interval  as  an  example,  the 
results  show  that,  in  Sample  1,  there  are  2  particles,  or  25%  of 
the  total  particles  that  have  perimeters  between  90  and  120 
microns,  and  that  50%  of  the  total  particles  have  perimeters  less 
than  120  microns.  The  data  for  the  4th  interval  for  Sample  2  show 
that  the  average  number  of  particles  in  the  two  images  with  per- 
imeters between  90  and  120  microns  is  1.5,  the  average  percentage 
of  particles  with  perimeters  between  90  and  120  microns  is 
17.05%,  and  the  average  percentage  of  particles  with  perimeters 
less  than  120  microns  is  43.18%.  Note  that  the  number  of  classes 
in   Sample  2  is  greater  than  that  in  Sample  1  because  the  biggest 
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perimeter  in  Sample  2  is  greater  than  that  in  Sample  1. 


Tables  3  to  5  show  the  distributions  of  maximum  chords  and 
areas.  Again,  the  number  frequency,  the  relative  frequency  and 
the  cumulative  frequency  are  included.  The  structure  of  these 
tables  is  similar  to  Table  2.  In  the  distribution  of  west-east 
maximum  chords,  the  interval  between  classes  is  50  microns.  In 
the  distribution  of  north-south  maximum  chords,  the  interval  is 
100  microns.  In  the  distribution  of  areas,  the  interval  is  A  x 
10   square  microns. 
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Table  1 

Sample  1 

The  length  unit  is  microns 

In  the  ()  is  the  mean  value  for  all  the  samples  above 


Image  size  (x   e8): 

Area  fraction  of  intersections: 

Total  Number  of  intersections: 
#  of  intersection/unit  area  (x 


e-8) 


Max.  perimeter  of  intersections  (x  e4) 
Min.  perimeter  of  intersections: 
Mean  perimeter  of  intersections: 

Maximun  area  of  intersections  (x  e8): 
Minimun  area  of  intersections: 
Mean  area  of  intersections: 

Max.w-e  extent  of  intersections  (x  e4) 
Min.w-e  extent  of  intersections: 
Mean  w-e  extent  of  intersections: 

Max.n-s  extent  of  intersections  (x  e4) 
Min.n-s  extent  of  intersections: 
Mean  n-s  extent  of  intersections: 


0.0355 
0.3057 
0.3057) 

8 
225.51 

225.51) 
0.1860 
0670 
1220 
1220) 
0026 
,0003 
,0014 
0014) 
,0590 
,0110 
,0371 
,0371) 
,0770 
,0180 
,0426 


(   0.0426) 


Labels 

Per ime  t e  r 

(x  e4) 

0 

0.0750 

2 

0.0670 

3 

0. 1160 

4 

0.  1240 

5 

0.  1860 

6 

0. 1600 

7 

0.1360 

8 

0.1120 

Area 
(x  e8) 
0.0003 
0004 
0010 
0016 
0026 
0026 
0014 
0010 


W-E  Extent 
(x  e4) 
0.0110 
0250 
0380 
0450 
0420 
0590 
0410 


0.0360 


N-S  Extent 
(x  e4) 
0.0330 
0.0180 
0.0320 
0.0480 
0.0770 
0.0570 
0.0420 
0.0340 
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Table  1  (continued) 

Sample  2 

The  length  unit  is  microns 

In  the  ()  is  the  mean  value  for  all  the  samples  above. 


Image  size  (x   e8): 

Area  fraction  of  intersections: 

Total  Number  of  intersections: 

#  of  intersection/unit  area  (x  e-8) 

Max.  perimeter  of  intersections  (x 
Min.  perimeter  of  intersections: 
Mean  perimeter  of  intersections: 

Maximun  area  of  intersections  (x  el 
Minimun  area  of  intersections: 
Mean  area  of  intersections: 

Max.w-e  extent  of  intersections  (x 
Min.w-e  extent  of  intersections: 
Mean  w-e  extent  of  intersections: 

Max.n-s  extent  of  intersections  (x 
Min.n-s  extent  of  intersections: 
Mean  n-s  extent  of  intersections: 


( 
e4)  : 


e4) 


e4) 


0.0355 
0.4798 
0.3927) 

11 
309.56 
267.54) 
0.2580 
0680 
1328 
1274) 
0041 
,0002 
0015 
0015) 
,0600 
,0120 
,0425 
,0398) 
,  1020 
,0130 
,0448 


(   0.0437) 


Labe Is 

Per ime  te  r 

(x  e4) 

0 

0.0730 

2 

0.  1240 

3 

0. 1300 

4 

0. 1140 

5 

0. 1200 

6 

0.2580 

7 

0.1910 

8 

0.  1410 

9 

0.1610 

10 

0.0680 

11 

0.0810 

Area 
(x  e8) 
0.0004 
0008 
0011 
0012 
0014 
0041 
0034 
0017 
0023 
0002 
0004 


W-E  Extent 
(x  e4) 
0.0180 
0.0600 
0.0510 
0.0360 
0.0470 
0.0510 
0.0580 
0.0410 
0.0540 
0.0120 
0.0400 


N-S  Extent 
(x  e4) 
0.0280 
0160 
0270 
0380 
0390 
1020 
0740 
0570 
0660 
0330 
0130 
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Table  2 
Distribution  of  Perimeters  of  Quartz 


The  length  unit  is  microns. 

In  the  ()  are  the  mean  values  for  all  the  samples  above. 


Sample 


Start-with 
(x  e4) 


Frequency 


Rel. Freq 


Cum. Freq . 


0.000 

0 

0.0000 

0.0000 

0)    ( 

0.0000)   ( 

0.0000) 

0.030 

0 

0.0000 

0.0000 

0)    ( 

0.0000)   ( 

0.0000) 

0.060 

2 

0.2500 

0.2500 

2)    ( 

0.2500)   ( 

0.2500) 

0.090 

2 

0.2500 

0.5000 

2)    ( 

0.2500)   ( 

0.5000) 

0.  120 

2 

0.2500 

0.7500 

2)    ( 

!   0.2500)   ( 

0.7500) 

0.  150 

1 

0. 1250 

0.8750 

:            i)      < 

0.1250)   ( 

;   0.8750) 

0.  180 

i 

0.  1250 

1.0000 

:            i)      < 

;   0.1250)   ( 

;   1.0000) 

2 

0.000 

0 

0.0000 

0.0000 

:   0.0000)    i 

:   0.0000)   1 

;   0.0000) 

2 

0.030 

0 

0.0000 

0.0000 

:   0.0000)    1 

;   0.0000)   < 

:   0.0000) 

2 

0.060 

3 

0.2727 

0.2727 

;   2.5000)    1 

;   0.2614) 

:   0.2614) 

2 

0.090 

1 

0.0909 

0.3636 

:   1.5000) 

;   0.1705) 

:   0.4318) 

2 

0.120 

4 

0.3636 

0.7273 

;   3.0000) 

[       0.3068) 

!   0.7386) 

2 

0.  150 

1 

0.0909 

0.8182 

:    i.oooo) 

{       0.1080) 

:   0.8466) 

2 

0.180 

1 

0.0909 

0.9091 

{       1.0000) 

[   0.1080) 

{       0.9545) 

2 

0.210 

0 

0.0000 

0.9091 

C   0.0000) 

(   0.0000) 

[       0.9545) 

2 

0.240 

1 

0.0909 

1.0000 

(   0.5000) 

(   0.0455) 

C   1.0000) 
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Table  3 
Distribution  of  Max.  W-E  Chords  of  Quartz 


The  length  unit  is  microns. 

In  the  ()  are  the  mean  values  for  all  the  samples  above. 


Sample 

t       Start-with 
(x  e4) 

Frequency 

Rel. Freq . 

Cum. Freq . 

0.000 

0 

0.0000 

0.0000 

0)    ( 

0.0000)   ( 

0.0000) 

0.005 

0 

0.0000 

0.0000 

0)    ( 

0.0000)   ( 

0.0000) 

0.010 

1 

0.  1250 

0. 1250 

1)    ( 

0.1250)   ( 

0. 1250) 

0.015 

0 

0.0000 

0.  1250 

0)    ( 

0.0000)   ( 

0. 1250) 

0.020 

0 

0.0000 

0. 1250 

0)    ( 

0.0000)   ( 

;   0.1250) 

0.025 

1 

0.  1250 

0.2500 

:            i)      < 

;   0.1250)   1 

;   0.2500) 

L     0.030 

0 

0.0000 

0.2500 

:            o)      i 

;  o.oooo) 

;   0.2500) 

L     0.035 

2 

0.2500 

0.5000 

:            2)      i 

;   0.2500) 

:   0.5000) 

L     0.040 

2 

0.2500 

0.7500 

:            2) 

C   0.2500) 

(   0.7500) 

L     0.045 

i 

0.  1250 

0.8750 

c            i) 

(   0.1250) 

(   0.8750) 

L     0.050 

0 

0.0000 

0.8750 

C       0) 

(   0.0000) 

(   0.8750) 

1     0.055 

1 

0.  1250 

1.0000 

C       i) 

(   0.1250) 

(   1.0000) 
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Table  3  (continued) 
Distribution  of  Max.  W-E  Chords  of  Quartz 


Sample  # 

Start-with 
(x  e4) 

Frequency 

Rel. Freq. 

Cum. Freq . 

2 

0.000 

0 

0.0000 

0.0000 

0.0000)    ( 

0.0000)   ( 

0.0000) 

2 

0.005 

0 

0.0000 

0.0000 

0.0000)    ( 

0.0000)   ( 

0.0000) 

2 

0.010 

1 

0.0909 

0.0909 

1.0000)    < 

0.1080)   ( 

0.1080) 

2 

0.015 

1 

0.0909 

0.  1818 

0.5000)    < 

;   0.0455)   ( 

0.  1534) 

2 

0.020 

0 

0.0000 

0.  1818 

;   0.0000)    ( 

:       0.0000)   ( 

;   0.1534) 

2 

0.025 

0 

0.0000 

0.  1818 

;   0.5000)    1 

;   0.0625)   < 

;   0.2159) 

2 

0.030 

0 

0.0000 

0.  1818 

;     o.oooo) 

;       0.0000) 

;   0.2159) 

2 

0.035 

1 

0.0909 

0.2727 

;   1.5000) 

:   0.1705) 

(   0.3864) 

2 

0.040 

2 

0.1818 

0.4545 

:   2.0000) 

(   0.2159) 

{       0.6023) 

2 

0.045 

1 

0.0909 

0.5455 

{        1.0000) 

(   0.1080) 

(   0.7102) 

2 

0.050 

3 

0.2727 

0.8182 

C   1.5000) 

(   0.1364) 

C   0.8466) 

2 

0.055 

1 

0.0909 

0.9091 

(   1.0000) 

(   0.1080) 

(   0.9545) 

2 

0.060 

1 

0.0909 

1.0000 

(   0.5000) 

(   0.0455) 

(   1.0000) 
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Table  4 
Distribution  of  Max.  N-S  Chords  of  Quartz 


The  length  unit  is  microns. 

In  the  ()  are  the  mean  values  for  all  the  samples  above. 


Sample 


Start-with 
(x  e4) 


Frequency 


Rel  .Freq 


Cum. Freq 


0.000 

0 

0.0000 

0.0000 

0)    ( 

0.0000)   ( 

0.0000) 

0.010 

1 

0.1250 

0.1250 

1)    ( 

0.1250)   ( 

0.1250) 

0.020 

0 

0.0000 

0.1250 

0)    ( 

0.0000)   ( 

0.1250) 

0.030 

3 

0.3750 

0.5000 

3)    ( 

0.3750)   ( 

0.5000) 

0.040 

2 

0.2500 

0.7500 

2)    ( 

0.2500)   ( 

0.7500) 

0.050 

1 

0.  1250 

0.8750 

1)    ( 

0.1250)   ( 

0.8750) 

0.060 

0 

0.0000 

0.8750 

0)    ( 

0.0000)   ( 

0.8750) 

0.070 

1 

0.  1250 

1.0000 

1)      < 

0.1250)   I 

1.0000) 

2 

0.000 

0 

0.0000 

0.0000 

0.0000)    ( 

;  o.oooo)  ( 

:   0.0000) 

2 

0.010 

2 

0.  1818 

0.  1818 

;   1.5000)    1 

;   0.1534) 

[   0.1534) 

2 

0.020 

2 

0. 1818 

0.3636 

;   1.0000) 

;   0.0909) 

(   0.2443) 

2 

0.030 

3 

0.2727 

0.6364 

:   3.0000) 

C   0.3239) 

(   0.5682) 

2 

0.040 

0 

0.0000 

0.6364 

(   1.0000) 

(   0.1250) 

(   0.6932) 

2 

0.050 

1 

0.0909 

0.7273 

(   1.0000) 

(   0.1080) 

(   0.8011) 

2 

0.060 

1 

0.0909 

0.8182 

(   0.5000) 

(   0.0455) 

(   0.8466) 

2 

0.070 

1 

0.0909 

0.9091 

(   1.0000) 

(   0.1080) 

(   0.9545) 

2 

0.080 

0 

0.0000 

0.9091 

(   0.0000) 

(   0.0000) 

(   0.9545) 

2 

0.090 

0 

0.0000 

0.9091 

(   0.0000) 

(   0.0000) 

(   0.9545) 

2 

0.  100 

1 

0.0909 

1.0000 

(   0.5000) 

(   0.0455) 

(   1.0000) 
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Table  5 
Distribution  of  Areas  of  Quartz 


The  length  unit  is  microns. 

In  the  ()  are  the  mean  values  for  all  the  samples  above. 


Sample  it       Start-with 
(x  e6) 


Frequency 


Rel. Freq 


Cum. Freq . 


0.000 

2 

0.2500 

0.2500 

2)    ( 

0.2500)   ( 

0.2500) 

0.040 

0 

0.0000 

0.2500 

0)    ( 

0.0000)   ( 

0.2500) 

0.080 

2 

0.2500 

0.5000 

2)    ( 

0.2500)   ( 

0.5000) 

0.  120 

2 

0.2500 

0.7500 

2)    ( 

0.2500)   ( 

0.7500) 

0.  160 

0 

0.0000 

0.7500 

0)    ( 

0.0000)   ( 

0.7500) 

0.200 

0 

0.0000 

0.7500 

0)    ( 

0.0000)   ( 

0.7500) 

0.240 

2 

0.2500 

1.0000 

2)    < 

0.2500)   < 

;   1.0000) 

2 

0.000 

3 

0.2727 

0.2727 

2.5000)    ( 

.   0.2614)   1 

;   0.2614) 

2 

0.040 

1 

0.0909 

0.3636 

;   0.5000)    < 

;   0.0455) 

;   0.3068) 

2 

0.080 

2 

0.  1818 

0.5455 

;  2.0000)   < 

;   0.2159) 

;   0.5227) 

2 

0.120 

1 

0.0909 

0.6364 

:   1.5000)    1 

;   0.1705) 

(   0.6932) 

2 

0.  160 

1 

0.0909 

0.7273 

;   0.5000) 

:   0.0455) 

{       0.7386) 

2 

0.200 

1 

0.0909 

0.8182 

(   0.5000) 

{       0.0455) 

(   0.7841) 

2 

0.240 

0 

0.0000 

0.8182 

(   1.0000) 

C   0.1250) 

(   0.9091) 

2 

0.280 

0 

0.0000 

0.8182 

(   0.0000) 

(   0.0000) 

(   0.9091) 

2 

0.320 

1 

0.0909 

0.9091 

(   0.5000) 

(   0.0455) 

(   0.9545) 

2 

0.360 

0 

0.0000 

0.9091 

(   0.0000) 

(   0.0000) 

(   0.9545) 

2 

0.400 

1 

0.0909 

1.0000 

(   0.5000) 

(   0.0455) 

(   1.0000) 
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Typical  Results  of  Estimation  of  3-Dimensional  Properties 

If  Figures  1  and  2  were  cut  planes  through  a  solid  sample, 
the  results  in  Tables  1  to  5  could  be  used  to  obtain  estimates  of 
the  properties  of  the  solid.  In  what  follows,  it  will  be  assumed 
that  Figures  1  and  2  are  such  planes.  (This  ruse  is  used  to  cir- 
cumvent the  need  for  preparing  a  special  sample  merely  to  demon- 
strate the  AQIA  package.)  The  3-dimensional  properties  that  can 
be  estimated  by  AQIA  are  shown  in  Tables  6  and  7.  These  two 
tables  also  show  the  length  unit  used  in  the  estimations  and  the 
mean  values  from  the  two  separate  images. 

Table  6  shows  the  volume  fraction,  the  surface  area  of 
objects  per  unit  bulk  volume,  and  the  surface  area  of  objects  per 
unit  object  volume  (based  on  the  former  two  values).  The  total 
number  of  objects  per  unit  bulk  volume  is  given  and  is  calculated 
from  the  particle  size  distribution.  Table  6  also  shows  the  mean 
diameter  of  the  objects,  assuming  that  the  objects  are  spheres. 


Table  7  shows  the  particle  size  distribution.  The  number 
frequency,  the  relative  frequency  and  the  cumulative  frequency 
are  included.  The  interval  between  classes  is  100  microns.  The 
mean  lengths  of  the  intercepted  chords  in  each  class  are  used  as 
the  mean  diameters  of  objects  in  that  class.  Table  7  also 
includes  the  number  frequency  of  intersected  chords'  lengths. 
For  example,  the  data  for  the  4th  interval  for  Sample  1  show  that 
there  are  109  intercepted  chords  with  mean  length  350  microns,  or 
8.6025  x  10~10  particles  with  a  mean  diameter  of  350  microns   per 
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Table  6 

Volume  Fraction  and  Surface  Areas 

VFA  is  volume  fraction.  < 

SV  is  surface  area  per  unit  bulk  volume. 

SSA  is  specific  surface  area  per  unit  object  volume. 

TN  is  the  total  number  of  particles  per  unit  bulk,  volume 

MD  is  the  mean  diameter  of  particles 

The  length  unit  is  microns 

In  the  ()  are  the  mean  values  for  all  the  samples  above. 


Sample  # 


1 


VFA 


SV(x  e-4)  SSA(x  e-4)  TN(x  e-10)  MD(x  e4) 


0.3057  29.9366  97.9280  44.8639    0.0367 

(0.3057)  (  29.9366)  (  97.9280)  (  44.8639)  (0.0367) 

0.4798  46.2093  96.3094  116.1319    0.0278 

(0.3927)  (  38.0729)  (  97.1187)  (  80.4979)  (0.0322) 


cubic  micron  of  bulk  volume,  the  percentage  of  particles  with  a 
mean  diameter  of  350  microns  is  19.17%,  and  the  percentage  of 
particles  with  diameter  less  than  400  microns  is  56.89%. 
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Table  7 


Particles  Size  Distribution 

The  length  unit  is  microns.  Chd  stands  for  intersected  chords. 
In  the  ()  are  the  maen  values  for  all  the  samples  above. 
Frequency  is  the  number  of  particles  per  unit  bulk  volume. 


Sample 


Ave . Dia . 
(x  e4) 


Freq.(x  e-10)   Rel.freq.    Cum.Freq.   Chd.Freq 


0.0050 

0.0000 

0.0000 

0.0000 

25 

0.0000)    ( 

0.0000)   ( 

0.0000) 

( 

25) 

0.0150 

14.7153 

0.3280 

0.3280 

75 

14.7153)    ( 

0.3280)   ( 

0.3280) 

( 

75) 

0.0250 

2.2047 

0.0491 

0.3771 

84 

2.2047)    ( 

0.0491)   ( 

0.3771) 

( 

84) 

0.0350 

8.6025 

0.  1917 

0.5689 

109 

8.6025)    ( 

0.1917)   < 

0.5689) 

( 

109) 

0.0450 

6.6163 

0.  1475 

0.7164 

97 

6.6163)    ( 

0.1475)   ( 

0.7164) 

( 

97) 

0.0550 

8.8598 

0.  1975 

0.9138 

78 

8.8598)    < 

0.1975)   ( 

0.9138) 

( 

78) 

0.0650 

0.9939 

0.0222 

0.9360 

28 

0.9939)    ( 

0.0222)   ( 

0.9360) 

( 

28) 

0.0750 

2.4490 

0.0546 

0.9906 

24 

2.4490)    < 

0.0546)   < 

;   0.9906) 

( 

24) 

0.0850 

0.4222 

0.0094 

1.0000 

4 

0.4222)    ( 

;   0.0094)   < 

;   1.0000) 

( 

4) 

2 

0.0050 

48.9697 

0.4217 

0.4217 

55 

.  24.4849)    1 

;   0.2108) 

;   0.2108) 

( 

40.00) 

2 

0.0150 

1. 1944 

0.0103 

0.4320 

83 

;   7.9548)    1 

;   0.1691)   t 

;   0.3800) 

( 

79.00) 

2 

0.0250 

0.0000 

0.0000 

0.4320 

135 

;   1.1023) 

:   0.0246)   t 

{       0.4045) 

(109.50) 

2 

0.0350 

41.5190 

0.3575 

0.7895 

247 

;  25.0608) 

{       0.2746) 

C   0.6792) 

(178.00) 

2 

0.0450 

0.3619 

0.0031 

0.7926 

109 

[   3.4891) 

(   0.0753) 

(   0.7545) 

(103.00) 

2 

0.0550 

19.2688 

0.1659 

0.9585 

131 

[     14.0643) 

C   0.1817) 

C   0.9362) 

(104.50) 

2 

0.0650 

0.0000 

0.0000 

0.9585 

15 

(   0.4970) 

(   0.0111) 

(   0.9473) 

( 

21.50) 

2 

0.0750 

2.3536 

0.0203 

0.9788 

25 

(   2.4013) 

(   0.0374) 

(   0.9847) 

( 

24.50) 

2 

0.0850 

0.  1609 

0.0014 

0.9802 

6 

(   0.2915) 

(   0.0054) 

(   0.9901) 

( 

5.00) 

2 

0.0950 

0.0000 

0.0000 

0.9802 

5 

(   0.0000) 

(   0.0000) 

(   0.9901) 

( 

2.50) 

2 

0.  1050 

2.3035 

0.0198 

1.0000 

27 

(   1.1517) 

(   0.0099) 

(   1.0000) 

( 

13.50) 
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CONCLUSIONS 

1.   An  automatic,  quantitative  image  analysis  (AQIA)  system   has 
been  developed  for  construction  materials. 


2.  A  special  technique  of  image  segmentation,  based  on  the  con- 
cepts of  fuzzy  probability,  has  been  developed.  It  models 
the  cognition  process  of  humans,  and  produces  consistent 
binary  images . 

3.  The  system  uses  object  labeling  to  process  the  binary  images 
more  accurately  than  non-labeling  image  analysis  systems. 

4.  In  the  case  of  discrete  particles,  the  system  can  measure 
the  perimeters,  areas,  maximum  chords  of  the  particles,  the 
maximum,  minimum  and  mean  values  of  these  parameters,  the 
distributions  of  these  parameters,  the  area  fraction  of  the 
particles  and  the  total  number  of  particles. 

5.  In  the  case  of  a  cut  surface  through  a  massive  sample,  the 
system  can  estimate  the  volume  fraction  of  objects,  the  sur- 
face area  of  objects  per  unit  bulk  volume,  and  the  surface 
area  of  objects  per  unit  object  volume.  Assuming  that  the 
objects  are  spheres,  the  particle  size  distribution  of 
objects  in  a  massive  sample,  the  total  number  of  objects  per 
unit  bulk  volume  and  the  mean  diameter  can  also  be 
e  s  t imat  ed . 


